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Abstract 

We consider spatial stochastic models, which can be applied e.g. to telecom- 
munication networks with two hierarchy levels. In particular, we consider Cox 
processes Xl and Xh concentrated on the edge set of a random tessel- 
lation T, where the points Xl 7 ti and Xh,u of Xl and Xh can describe the 
locations of low-level and high-level network components, respectively, and 
the underlying infrastructure of the network, like road systems, railways, 
etc. Furthermore, each point X^^ of Xl is marked with the shortest path along 
the edges of T to the nearest (in the Euclidean sense) point of Xh- We investi- 
gate the typical shortest path length C* of the resulting marked point process, 
which is an important characteristic e.g. in performance analysis and planning 
of telecommunication networks. In particular, we show that the distribution 
of C* converges to simple parametric limit distributions if a scaling factor k 
converges to zero and infinity, respectively. This can be used to approximate 
the density of C* by analytical formulae for a wide range of k. 
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1 Introduction 

Asymptotic properties of spatial stochastic models are considered, which can be ap- 
plied e.g. in the analysis and planning of telecommunication networks. More precisely, 
we consider stochastic models for networks with two hierarchy levels, i.e., there are 
network components of two different kinds: low-level components (LLC) and high- 
level components (HLC). The locations of both HLC and LLC are represented by 
points in the Euclidean plane M. 2 . We then associate with each HLC a certain subset 
of R 2 which is called its serving zone. This is done in such a way that the serving 
zones of the HLC are disjoint convex polygons which cover the whole R 2 . Each LLC is 
linked to the HLC in whose serving zone the LLC is located. In particular, we assume 
that the serving zones are constructed as the cells of the Voronoi tessellation with 
respect to the locations of HLC. This is equivalent to link each LLC to its nearest 
HLC, where "nearest" means with respect to the Euclidean distance. Furthermore, 
we assume that the HLC and LLC are located on the edges of a random geometric 
graph, where the link from a LLC to its nearest HLC is assumed to be the shortest 
path along the edges of that graph. In the case of telecommunication networks the 
edges of the random geometric graph represent the underlying infrastructure, e.g. an 
inner-city street system. 

Thus, we study a class of stochastic network models which has been introduced 
in [10] as the Stochastic Subscriber Line Model (SSLM) for urban access networks. 
Note that the SSLM is a model from stochastic geometry which provides tools for 
the description of geometric features of the network. Based on this model, stochastic 
econometrical analysis can be done for real telecommunication networks, e.g. connec- 
tion costs for access networks can be determined, see [EUdU 36, 38J, where we focus 
on the case that the infrastructure of the network is modeled by the edge set of a 
stationary random tessellation and both the HLC and LLC are modeled by Cox pro- 
cesses concentrated on this edge set. Then we are especially interested in the shortest 
path length along the edge set between LLC and HLC, which is an important perfor- 
mance characteristic in cost and risk analysis as well as in strategic planning of wired 
telecommunication. In order to define an appropriately chosen (global) distribution 
of the shortest path length we regard the so-called typical shortest path length C* . 
It can be interpreted as the length of the shortest path from a location of LLC, which 
is chosen at random among all locations of LLC, and its nearest HLC. We are then 
interested in the asymptotic behaviour of the distribution of C* for two extreme cases 
of model parameters. In particular, we show that the distribution of C* converges 
to simple parametric limit distributions if a scaling factor k converges to zero and 
infinity, respectively. This can be used to approximate the density of C* by analytical 
formulae for a wide range of k which is a great advantage e.g. for the econometrical 
analysis of real telecommunication networks, see [H]. The mathematical techniques, 
which we exploit in order to derive our main results presented in Theorems 13.11 and 
13.21 include Palm calculus and Poisson approximation for stationary point processes, 
Kingman's subadditive ergodic theorem, and the generalized Blaschke-Petkantschin 
formula from geometric measure theory. 

The paper is organized as follows. In Section [2] we give a short description of 
the particular stochastic network model considered in the present paper. Then, in 
Section EH we present the main results stated in Theorems 13.11 and I3.2L The proof of 
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Theorem 13 .21 is given in Section Q], where some details are postponed to the Appendix. 
In Section [5], it is shown that the mixing and integrability conditions of Theorems 13. II 
and 13.21 are fulfilled for various examples of random tessellations. Some extensions 
of our results to other performance characteristics, more general classes of random 
geometric graphs, and more general connection rules are discussed in Section 15 .41 
Finally, Section[6]concludes the paper and gives an outlook to possible future research. 

2 Stochastic modelling of hierarchical networks 

To begin with we give a short description of the particular stochastic network model 
considered in the present paper. For more details on this model see also [T3]. More- 
over, we briefly explain the mathematical background and introduce the notation we 
are using. For further details on spatial point processes and random tessellations, 
see e.g. [HI EH [32], 133] . Surveys on applications of tools from stochastic geometry 
to spatial stochastic modelling of telecommunication networks can be found e.g. in 

[HISS]. 

2.1 Marked point processes 

First we recall some basic notions and results regarding marked point processes in M 2 . 
They can be used to model locations of customers or equipments in telecommunication 
networks. Let B 2 denote the family of Borel sets of M 2 and N the family of all simple 
and locally finite counting measures on B 2 . Note that each <p G N can be represented 
by the sequence {x n } of its atoms, i.e. <p = ^2 n S Xn , where S x is the Dirac measure 
with S X (B) — 1 if x G B and S X (B) = if x B. Let Af denote the cx-algebra of 
subsets of N generated by the sets {ip £ N : p(B) = j} for j G N and B G B 2 . The 
shift operator t x : N i-> N is defined by t x p{B) = <p(B + x) for x G M 2 and B G B 2 , 
where B + x = {x + y : y <E B}. Then a point process X is a random element of 
the measurable space (N,J\f), where we identify X with the sequence {X n } of its 
(random) atoms, writing X = {X n } for brevity. 

Let M be a Polish space with its Borel cx-algebra Bm- Then we use the notation 
N M for the family of all counting measures on B 2 ® Bm which are simple and locally 
finite in the first component. Note that the atoms (x n ,m n ) of the counting measure 
if) = ^2 n 8(x n ,m n ) £ -Wm have two components: the location x n G M 2 and the mark 
m n G M. The cx-algebra A/m is defined in the same way as above and the shift 
operator t x : N M \— > N M translates the first component of the atoms of ip G iV M by 
-x, i.e. t x (if)) = J2n S (xn-x,mn)- A random element X = {(X n , M n )} of (N M ,M M ) is 
then called a marked point process. 

2.2 Palm distributions 

Stationarity and isotropy of (marked) point processes are defined in the usual way, i.e., 
assuming the invariance of their distributions with respect to arbitrary translations 
and rotations around the origin, respectively. By A > we denote the intensity of 
a stationary marked point process X = {(X n , M n )}, i.e. A = E#{n : X n G [0, l] 2 }, 
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and the Palm mark distribution ¥° x : £>m — » [0, 1] of X is given by 

n(G) = E #{, i :X„ £ [0,l)'.M. £G}< ^ (21) 

A 

A random variable M* distributed according to F x is called the typical mark of X. 

Furthermore, two jointly stationary marked point processes X* 1 ) = {(Xn , Mn^)} 
and = {(Xn \ Mn^)} with intensities Ai and A 2 and mark spaces Mi and M 2 , 
respectively, will be considered as random element Y = (X^\ X^) of the product 
space X Ml M2 = N Ml x X M2 . The Palm distribution P^ w of Y with respect to the 
i-th component, i = 1, 2, is then defined on Ami ® Am 2 ® ^M, by 

E#{n : Xj? G [0, 1) 2 ,M^ G G,t x{i) Y G A} 
F* xii) (AxG) = — *= L j (2 . 2) 

where A G Ami ® Am 2 and G G £>M; ■ Note that the Palm mark distribution P^) of 
I™ can be obtained from P^j as a marginal distribution. 



2.3 Random tessellations 

As a model for the underlying random geometric graph we consider the edge set 
of random tessellations of M. 2 . Note that a random tessellation T is a locally finite 
partition {S n } of M 2 into random (compact and convex) polygons H n , which are called 
the cells of T. We can also regard T as a marked point process {(a(5 n ), 5°)}, where 
the shifted cells S° = S n — a(E n ) contain the origin. The points a(E n ) G S n C M 2 are 
then called the nuclei of the cells 5 n of T. Furthermore, we can identify T with the 
edge set T"W = [J n d'E n of T. Note that T^ 1 ' is a random closed set in M 2 , i.e., T^ 1 ' 
is a random element of (J 7 , £>(.F)), where T denotes the family of all closed subsets 
of M 2 and B{T} is the smallest cr-algebra of subsets of T which contains the „hitting 
sets" Tc — {B G T : B n C ± 0} for all compact C G B 2 . 

If T is stationary, i.e., T^> = T^ 1 ) + x for each x G M 2 , then the intensity 7 of 
T is defined as 7 = Ez/^T^ 1 ) fl [0, l] 2 ), i.e. the mean length of per unit area, 
where i>\ denotes the 1-dimensional Hausdorff measure. In the following we always 
assume that T is a (normalized) stationary tessellation with Kui(T^ fl [0, l] 2 ) = 1. 
Furthermore, for each 7 > we consider the scaled tessellation T 7 with intensity 7 
which is defined by T 7 = T/7, i.e., we scale the edge set with I/7 getting Tj 
such that Ez/i(T^ 1} n [0, l] 2 ) = 7. 

A random tessellation T is called isotropic if the distribution of is invariant 
with respect to rotations around the origin. Furthermore, a stationary tessellation T 
is called mixing if 

lim P(T (1) G A, T (1) + x G A') = P(T (1) G A) P(T (1) G A') 

\x\ >oo 

for any A, A' G B^J 7 ). Note that for any T which is mixing it holds that 

P(T (1) GA) = 1 or P(T (1) eA)=0 for each A G T{T) , (2.3) 

where X(JF) denotes the sub-a-algebra of invariant sets of B(F), i.e. A + x = A for 
all A G X(jF) and x G IR 2 . A stationary tessellation T which satisfies condition (12.31) 
is said to be ergodic. 
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2.4 Cox processes on edge sets 

For any 7 > 0, we consider Cox point processes Xh = {Xn,n} and Xl = {^L,n} 
concentrated on T 7 , in order to model the locations of HLC and LLC, respectively. 
In particular, we assume that Xh is a Cox process on T 7 with linear intensity \i 
which is constructed by placing homogeneous Poisson processes on the edges of T 7 
with linear intensity Xg. The random driving measure Ax H '■ B 2 — > [0, 00] of Xh is 
then given by 

A IjJ (5) = V 1 (BnTW), BeB 2 . (2.4) 

Analogously, X L is a Cox process on with linear intensity X' e which is constructed 
in the same way, i.e., by placing Poisson processes on the edges of T 7 with linear 
intensity \' t Thus, Xh and Xl are Cox processes concentrated on the same edge 
set where we assume that X H and X L are conditionally independent given T 7 . 
Furthermore, note that X H and X L are stationary, isotropic, and ergodic if T is 
stationary, isotropic, and ergodic, respectively. The planar intensities A and A' of Xh 
and Xl are given by A = A^7 and A' = A^y. 

2.5 Serving zones and shortest paths 

Let Th = {Sff>} denote the Voronoi tessellation induced by the points Xn, n °f the 
Cox process Xh = {Xn, n }, i-e. 

^H,n = {x E M 2 : \x — Xn,n\ < \x — XH, m \ f° r all m ^ n} , 

where | • | denotes the Euclidean norm. The cells of Th are considered to be 
the serving zones of HLC. By means of the four modelling components T 7 , Xh, Xl 
and Th we can construct the marked point process Xl,c — {(XL,n>C n )}, where the 
mark C n is the length of the shortest path from Xl, u to Xhj along the edge set T 7 ^ 
of T 7 provided that Xl,u £ ^hj- 

Thus, each LLC is connected to its nearest HLC in the Euclidean sense and not 
in the shortest path sense. However, for applications this is a reasonable assump- 
tion since the planning of telecommunication networks is complicated and existing 
networks have evolved for long periods. Therefore, it is unrealistic to assume that 
serving zones are defined with respect to the shortest path distance and it is appro- 
priate to use a simpler rule. Furthermore, analysis of real data has shown that the 
approach considered in the present paper is realistic ([14]). 

It is not difficult to show that X L) c is a stationary and isotropic marked point 
process if T 7 is stationary and isotropic, respectively. Realizations of service zones 
and shortest paths are displayed in Figure [TJ(a) and (b) for T 7 being a Poisson- Voronoi 
tessellation (PVT) and a Poisson line tessellation (PLT), respectively. 

The model characteristic we are mainly interested in is the distribution of the 
typical mark C* of Xl@- Thus, we are interested in the Palm mark distribution 
^"x LC °f Xl,Ci i- e -i the distribution of the typical shortest path length. 

Note that the realizations of Xl,c can be constructed from the corresponding 
realizations of Xl and Xh,s, where Xh,s = {(^H,n, S H n )} is a stationary marked 
point process with marks S H n = (T^ nS^„) —X H , n - Thus, instead of X Lj c, we can 
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(a) PVT as infrastructure model (b) PLT as infrastructure model 

Figure 1: Higher-level components with their serving zones (black) and lower-level 
components (grey with black boundary) with shortest paths (dashed) along the edge 
set (grey). 

consider the vector Y = (Xl, Xh,s) and the Palm distribution of Y with respect 
to Xl, which has been introduced in (12.21) . Let (X^,Xh,s) be distributed according 
to P^, where we use the notation X H S = {pCe>, S 1 ^ n )} and 

T« = |J (S° H>n + X H , n ) . (2.5) 

n>l 

Note that X H = {X H n } is a (non-stationary) Cox process on Tj with linear intensity 
A^. Moreover, by Xh,o we denote the closest point (in the Euclidean sense) of {XH, n } 
to the origin. Then, the typical shortest path length C* can be given by C* = c(X H fl), 
where c(X Hj0 ) denotes the length of shortest path from the origin to X Hj0 , along the 
edges^of In the following we always assume that the joint distribution of C*, Xh 
and T 7 is given by F* Xl . 

3 Limit theorems for the typical shortest path length 

We investigate the asymptotic behavior of the distribution of C* for two different 
cases: 7 — > with fixed and •y/Xi — > 00 with •yXi fixed, i.e., unboundedly sparse 
edge sets and unboundedly dense edge sets, respectively. For 7 — » 0, we show in The- 
orem [3J] that the distribution of C* converges weakly to an exponential distribution, 
where no specific assumption on the underlying stationary tessellation T is needed. 
Furthermore, for 7 — > 00 and T being a stationary and isotropic random tessellation 
which is mixing, we get in Theorem 13.21 that the distribution of C* converges weakly 
to a Weibull distribution. 
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(a) k = 0.5 (b) k = 1000 

Figure 2: Realizations of Xh,s — {(^H,n, S# n )} for extremal values of k 

3.1 Scaling invariance property 

Recall that the stochastic network model introduced in Section [2] and, in particular, 
the distribution of C* is fully specified by T, 7, A^ and X' e . Moreover, it can be shown 
(see e.g. [EH [36]) that the distribution of C* does not depend on X' e . Therefore, we 
only regard the parameters 7 and Xi in the following. Sometimes we use the notation 
C* = (7* (7, A|) to emphasize that the distribution of C* depends on 7 and A^. 

Furthermore, a scaling invariance property holds for this model. If the value of the 
quotient K = 7/A^ is constant, then the structure of Xh,s is fixed, but on different 
scales for different parameter vectors (7, A^) = Xg). We are interested in the 
limiting behavior of the distribution of C* for k — > with A^ fixed and for k — > 00 
with A = A^7 fixed. In Figure [2] realizations of Xh,s are shown for two (extremely 
small and large) values of k, where the realization of T is sampled from a PLT. One 
can see that for small values of n the segment systems within the serving zones mainly 
consist of one single segment only, whereas for large values of k the networks inside 
the serving zones become rather dense. 



3.2 Asymptotic exponential distribution for k — > 

First we regard the case that k = 7/A^ — > with A^ fixed, i.e., 7^0. 

Theorem 3.1 Let T be an arbitrary stationary tessellation. Then, for any fixed 
Xe > 0, it holds that 



C*( 7 ,A,) 



as 7 



0. 



(3.1) 



where — > denotes convergence in distribution and Z ~ Exp(2A^), i.e., the random 
variable Z is exponentially distributed with expectation (2A^) _1 . 



Proof Let i? 7 = max{r > : B(o, r) fl L° = B(o, r) R T^}, where B(o, r) denotes 
the ball centered at the origin with radius r and L° is the segment containing the 
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origin of the random edge set introduced in (12.51) . It is not difficult to see that 

limi? 7 = oo a.s. (3.2) 

7— >o 

Recall that C* = c(X Hi0 ), where X H>0 is the closest point to the origin of the point 
process X H = {X Hn } of HLC under Px L , and note that the values of the distribution 
function F c * : (0, 00) — > (0, 1) of C* can be written as 

F c * (x) = F(X Hfi £ B(o, R 7 )) ¥(C* < x \ X Hfl G B(o, Ry)) 

+ F(X Hfi £ B(o, Ry)) F(C* < x I X Hfi £ B(o, Ry)) 



for each x > 0. It can be shown (see e.g. [9]) that X H is a Cox process which is 
a homogenec 
implies that 



a homogeneous Poisson process with linear intensity A^ on the edges of . This 



P(C* < x I X Hfi e B(o,Ry)) 



P(min{Zi, Z 2 } < x,min{Zi, Z 2 } < R 



1) 



F(X H ,o eB(o,R,)) 



for each x > 0, where the random variables Z-y and Z 2 are independent, exponentially 
distributed with parameter Xe and independent of i? 7 . Furthermore, we get that 

F(X H ,o g B(o, Ry)) = P(min{Z 1 , Z 2 } > Ry) = Eexp(-2A £j R 7 ) , 

since mm{Zx, Z 2 } is exponentially distributed with parameter 2A^ and independent 
of R-y. Thus, using (|3.2p . it follows that 

lim F(X Hfi & B{o, iU) = and lim F{X Hi0 € B{o, Ry)) = 1 
7— >o ' 7— >o ' 

and, consequently, lim 7 ^ F c * (x) = P(min{Z 1 , Z 2 } < x) = 1 — exp(— 2A^x) for each 
x > 0. 

Note that the case k = •y/Xe — > with 7 fixed and A^ — > 00 can be treated in the 
following way. Due to the scaling invariance property mentioned in Section 13.11 we 
have 

A, C*( 7 ,A,) = C*( 7 /A,,1) 
for any 7, A<? > 0. Thus, Theorem 13.11 yields that 

A^ C*(j, Xi) — > Z as Xe — > 00 , 

where Z ~ Exp(2). 

3.3 Asymptotic Weibull distribution for k — > 00 

In this section we assume that T is a stationary and isotropic random tessellation 
which is mixing. Furthermore, we assume that 

Euf(dE*) < 00, (3.3) 

where vi(&E*) denotes the circumference of the typical cell S* of T. 
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We investigate the asymptotic behavior of the distribution of C* = C*(7, Xg) for 
k — > oo, where 7 — ► 00 and A^ — ► such that A^7 = A is fixed. In particular, we show 
that C* converges in distribution to £Z, where £ > 1 is a certain constant which is 
multiplied by the (random) Euclidean distance Z from the origin to the nearest point 
of a stationary Poisson process of intensity A. Then, it is easy to see that Z as well 
as £Z have Weibull distributions. 

Theorem 3.2 Let Z ~ Wei(A7r,2) for some A > 0. Then there exists a constant 
£ > 1 such that 

C*( 7 ,A £ )^£Z as k^oo (3.4) 
provided that 7 — > 00 and Xi — > A^7 = A, where £Z ~ Wei(A7r/£ 2 , 2). 

The proof of Theorem 13.21 is split into several steps. We first show in Lemma l4~2l 
that under the Palm probability measure P^i the Euclidean distance |-Xfr,o| from the 
origin to the nearest point Xh,q of the point process X H = {X HtTl } of HLC converges 
in distribution to the corresponding characteristic of a stationary Poisson process 
with intensity A. Furthermore, in Lemma [4.41 we show that for some constant £ > 1 
the difference between £|X# i0 | and the shortest path length C* = C*(j, Xe) from 
the origin to X Hj0 along the edge set Ty converges in probability to zero. Then, 
combining the results of Lemmas 14.21 and I4T4"1 the assertion of Theorem 13.21 follows. 



4 Proof of Theorem 3.2 



4.1 Some auxiliary results on convergence of point processes 

In the proofs of Lemmas 14.11 and 14.21 which will be given below, we use two classic 
results regarding the convergence in distribution of point processes, see e.g. [HI IB, 26j. 
Note that a sequence of point processes X^\X^ 2 \ ... in R 2 is said to converge in 
distribution to a point process X in R 2 if 

lim P(X^( J B 1 ) =*!,.. . ,X^ m \B k ) = i k ) = F(X(B 1 ) = i x , . . .,X(B k ) = i k ) 

for any k > 1, ii, . . . , i k > and for all finite sequences of bounded sets B\, . . . , B k € 
B 2 which satisfy the condition ¥(X(dBj) > 0) = for each j — 1, . . . , k. In this case 
we shortly write X^ => X. 

Let X = {X n } be an arbitrary ergodic point process in R 2 with P(X(R 2 ) = 0) = 0, 
and let A G (0, oo) denote the intensity of X. Then, the following limit theorem for 
independently thinned and appropriately re-scaled versions of X is true. For each c G 
(0, 1), let X^ denote a point process which arises from X by independent thinning, 
where each atom X n of X is deleted with probability 1 — c (and „survives" with 
probability c). Furthermore, let Y^ c > be a re-scaled version of X^ c \ where Y^ C \B) = 
X^(B/y/E) for each B G B 2 . Then, for each c G (0,1), the point process is 
stationary with the same intensity A as X, and 

F( C )=^Y ifc^O, (4.1) 

where Y is a stationary Poisson process in R 2 with intensity A, see e.g. Section 11.3 
of [8j or Theorem 7.3.1 in [26]. Moreover, the following continuity property of Palm 
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distributions holds. Let X,X^' ,X^ 2 \ . . . be stationary point processes in IR 2 such 
that P(X(M 2 ) = 0) = P(X( m )(M 2 ) = 0) = for each m > 1 and let A, A (1) , A (2) , . . . 
denote the intensity of X, X^,X^ 2 \ . . ., respectively. If A^ m ) = A for each m > 1 and 
X( m ) ==^ x as m — > oo, then 

yM ^ Y as m -> oo, (4.2) 

where Y, Y Y ^ , . . . are point processes in M 2 whose distribution is equal to the 
Palm distribution of X, X^\ X^ 2 \ . . ., respectively, see e.g. Proposition 10.3.6 in [26]. 

4.2 Euclidean distance from the typical LLC to its closest 
HLC 

Throughout this section we assume that thejmderlying tessellation T is ergodic. In 
order to prove that the Euclidean distance |^Cff,o| from the typical LLC to its closest 
HLC is asymptotically Weibull distributed, we first show that the (stationary) Cox 
process X H converges in distribution to a homogeneous Poisson process if k — > oo 
provided that A^7 = A is constant. 

Lemma 4.1 If k —> oo, where \f) = A for some constant A G (0, oo), then Xh =^ 
Y , where Y is a stationary Poisson process with intensity A. 

Proof For each 7 > 1, let Xh = Xh(j) denote the Cox process of HLC with 
parameters 7 and A^, where A^ = A/7 for some constant A G (0, 00). Note that 
the Cox process X H {pf) can be obtained from X H {1) by independent thinning with 
survival probability c = I/7 and by subsequent re-scaling with scaling factor a/I/V 
Furthermore, the Cox process XniX) is ergodic, since T is ergodic. Thus, using (I4.1B . 
we get that Xh(j) =>• Y as 7 — > 00. 

Lemma 4.2 Let Z ~ Wei(A7r,2) for some A > 0. Then \Xh,q\ —> Z as k — > 00 
provided that 7 — > 00 and A^ — > such that = A. 

Proof Let X H {p/) be a point process in M 2 whose distribution is equal to the Palm 
distribution of X H = X H {pf). Furthermore, let Y be a stationary Poisson process 
with intensity A. Note that the distribution of Y + 5 Q is then equal to the Palm 
distribution of Y, see e.g. Proposition 13.1. VII in [8]. Thus, using (14.21) . Lemma I4~T1 
gives that 

X^^Y + So (4.3) 

as 7 — > 00 and A<? — > 0, where A^7 = A. Since X L and X H are Cox processes 
concentrated on which are conditionally independent given Tj , we get that 
Xh + d~o and the Palm version X H of Xh have the same distributions. This is an easy 
consequence of the representation formula for the Palm distribution of stationary Cox 
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processes, see e.g. Section 5.2 in [33]. In particular, this gives that for each r > 
lim F(\X H . \ > r) = lim F(X H (B(o, r)) = 0) 

= UmF((X H + 5 )(B(o,r)) = l) 

7^00 

= \im¥(X* H (B(o,r)) = 1) 

7— »oo 

= F((Y + 6 )(B(o,r)) = l) 
= P(l%B( o ,r))=0), 

where we used (14.31) in the last but one equality. Thus, for each r > 0, 
lim P(\Xh,o\ > r) = F(Y(B(o,r)) = 0) = exp(-A7rr 2 ) , 

7— >oo 

which means that |-X# )0 | — > ^ ~ Wei(A7r,2). 

4.3 Shortest path length vs. scaled Euclidean distance 

In this section we assume that T is a stationary and isotropic random tessellation 
which is mixing. Furthermore, we assume that the integrability condition (13.31) is 
satisfied. Then, we can show that for some constant £ > 1 the difference between 
£|Xff,o| an d the shortest path length C* = C*(7,A^) from the origin to X H0 along 
the edge set T^p converges in probability to zero. In order to show this we need the 
following auxiliary result. 

Lemma 4.3 LetTjty ={uG : |c(w)— £|-u|| < e}, where £ > 1 is some constant 
and c(u) denotes the length of the shortest path from u to the origin along the edges 
0/T7 . 1/7 — > 00 and Xf — > 0, where A^7 = A is fixed, then there exists £ > 1 such 
that for each e > and r > 

lim Eexp ( - - ^(f W\T« £ n 5(o,r))) = 1 • (4.4) 

The ^roo/ of this lemma is postponed to the Appendix. Now, using Lemma 14.31 
we are able to complete the proof of Theorem 13.21 by showing that the following is 
true. 

Lemma 4.4 If ^ — > 00 and Xe — > such that Xej = X, then there is a constant £ > 1 

~ p p 
with C*yy, Xe) — £\Xh,o\ — > 0, where — > denotes convergence in probability. 

Proof We have to show that there exists a constant £ > 1 such that for any e > 
and 5 > we can choose 70 > with 

P(|C*-£|X H> o|| >e) <5 

for all 7 > 7o- Note that 
F(\C*-t\X H>0 \\>e) 

= F(\C*-Z\X H ,o\\ >£,\Xh,o\ <r) +F(\C*- £\X H ,o\\ > e,\X H ,o\> r) , 

11 
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where r > is an arbitrary fixed number. Since 

f(\X Hfi \ >r) — ► e~ Xnr2 as 7 — > 00 , 

see Lemma [Ol we can choose r > such that P(|X^ | > r) < 5/2 for all 7 > 
sufficientlyjarge. Thus, it is enough to show that there exists 70 > such that 
P(|C*- Z\X H)0 \\ > e, \X Hfi \ < r) < 5/2 for all 7 > 70. Let N = X H (B(o,r)) denote 
the number of points of Xh in B(o, r). Then we have 

W(\C*- (i\X Hfl \\ >e,\X Hfi \<r) 

too 
^>(iV = k I f 7 ) P(.max fc (\c{Yi) - £\Yi\\) > e 
k=l % ~ 

too 
^P(iV = A;|f 7 ) (l-F(|c(H)-emi| <e | ^r) 
fc=i 

where the points Yi, . . .,Yk are conditionally independent and identically distributed 
according to z/i( ■ C\T^ n £>(o, r))/vi (T 7 n £>(o, r)) for given T 7 and N = k. In 
particular, for the conditional probability in the latter expression, we have 

l[- £ , e ](c(u) - f|u|) i/i(du) 



T y ,N=k 



P(|c(ll) - < e I 3V) 



^(rfnfl^r)) 

^ff, E nB( ,r)) 



d(T 7 (1) n%r)) 

Using that iV ~ Poi(X) with A = A^i(T 7 Pi 5(o, r)) given T 7 , we get 

00 

£>(j\r = fc|f 7 ) (i-p(| c (yi)-eN| <^ I 



fc=i 

x _ T A fe /A^i(T (1) 
e 

fc=0 

00 



^ ,]iAV ^(r^nfl( ,r)) ^ 
1 k\\ \ ) 



k=0 

_ l _ e -Af( i ,i(f^ ) nB( ,r))- i /i(f^^ £ nB( ,r))) 

Thus we have 



\imF(\C*-£\X Hfi \\ >e,\X H ,o\ <r) < 1- lim Eexp (-- uM^T^.nB^r)) 

7— >oo 7— >oo V 'V " s ' 



Using Lemma PI this gives that lim^oo P(|C* — £|X# >0 || > £ 3 |^h,o| < r ) = 0, 
which completes the proof. 
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5 Examples 

Recall that in Theorem l3.2l we assumed that the underlying tessellation T is stationary 
and isotropic. The examples of tessellations discussed in the present section obviously 
possess these properties. Furthermore, we assumed in Theorem 13.21 that T is mixing 
and fulfills the integrability condition (13.31) . We first show that the mixing condition 
is satisfied for a wide class of tessellations. Moreover, we also show that (13.31) is true 
for these tessellations. 

The tessellation models considered in the literature focus mainly on PLT and 
PVT as well as on Poisson-Delaunay tessellations (PDT), on iterated tessellations 
constructed from these basic tessellations of Poisson type and on STIT tessellations, 
see e.g. [1] [5j , [9] [H], [30] and [36]-[40]. Here, we assume that an iterated tessella- 
tion is either a Tj/T/j-superposition or a Tj/T/j-nesting of tessellations Tj and Tn as 
defined e.g. in [21 [25], [40]. Note that the edge set of a T//T/7-superposition is given 
by the union Tj 1 ' U Tjj , where Tj and Tn are independent. Furthermore, a Tj/Tjj 
nesting is constructed by subdividing each cell of Tj by independent copies of T u . 
We show that for these important models Theorem 13.21 can be applied. Furthermore, 
if T is a PLT or a Tf/T/7-superposition/nesting with Tj being a PLT, then we can 
even calculate the constant £ explicitly that appears in Theorem I3.2L On the other 
hand, if T is a PDT, we get an upper bound for £. 

5.1 Mixing tessellations 

In order to apply Theorem 13.21 we have to show that the underlying tessellation T is 
mixing, where we will use the following criterion to show that a stationary random 
closed set is mixing. 

Lemma 5.1 A stationary random closed set 5 in M 2 is mixing if and only if 

lim p(s n d = 0, s n (c 2 + x) = 0) = p(s n c x = 0) p(s n c 2 = 0) (5.1) 

|x|^oo 

for all Ci, C 2 G 7Z, where 7Z is the family of all subsets ofM? which are finite unions 
of closed balls with rational radii and centres with rational coordinates. 

Note that the statement of Lemma I5TT1 is essentially Lemma 4 in [17], see also 
Theorem 9.3.2 in [32], where the (stronger) condition is considered that (15. ID holds 
for all compact sets Ci, C 2 C M 2 . However, it is easy to see that it suffices to assume 
that (I5.ip holds for the separating class 71; see also Section 1.4 of [29]. To make 
this clear, we only have to show that S = {J 7 ^ Ck '■ C , . . . ,Cf. G 1Z', k > 0} is a 
semi-algebra which generates B{J-), where TZ' = 1Z U and 

^8l..,c k = {Fe ^:FnC o = 0,FnC 1 ^0,...,FnC fc ^0}. 

Note that the family 7Z' is union-stable. Thus, by Lemma 2.2.2 in [32], we get that 
£ is a semi- algebra. Moreover, let G C R 2 denote an open set, then G = U^i ^ 
for some C u C 2 , . . . G 7Z' and T G = {F G T : F n G ^ 0} = U^i-^UU^' thus 
Tq G o{£). Since {JF G :GcK 2 open} generates &{T\ we get that &(£) = £>(JF). 
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Now the statement of Lemma I5TT1 can be proven by exactly the same arguments used 
in the proof of Lemma 4 in |17| . 

It is well known that T is mixing if T is a PDT, PVT and PLT, respectively, 
see e.g. Chapter 10.5 in [32], and recently it was shown that STIT tessellations are 
mixing (|22j). Furthermore, using Lemma [5.11 we can show that T is mixing if T is 
an iterated tessellation constructed from these basic tessellations of Poisson type. 

Lemma 5.2 The tessellationT is mixing ifT is a Tj/Tjj -superposition of two mixing 
tessellations Tj and Tu, or a Ti /Tu -nesting of a mixing initial tessellation Ti and 
any stationary component tessellation Tu. 

Proof Suppose first that T is a T//T//-superposition. Then, for any d,d G 1Z 

P(T (1) n Ci = 0, T (1) n (d + x) = 0) 

= P(T, (1) n d = 0, n (d + x ) = 0, Tff n d = 0, r£ J n (C 2 + x) = 0) 
= p(t 7 (1) n d = 0, t 7 (1) n (d + x) = 0) P(T / ( / 1) n d = 0, r# n (d + x) = 0) , 

since Tj and Tjj are independent. Thus, using Lemma [5.11 we get that T is mixing 
if Ti and Tu are mixing. Let now T be a Tj/Tu- nesting and assume that C\ = 
Uj =1 Bj,C 2 = lJ";=™ +1 Bj for closed balls B x , . . . , B n+m C M 2 with rational radii and 
centres with rational coordinates. Let S 1; S 2 , . . . be the cells of the initial tessellation 
Ti = {S n }, let D denote the family of all decompositions of the index set {1, . . . , n + 
m} into nonempty subsets, and for J = { Ji, . . . , J^} G D consider the set 

Aj(x) = {U jeJi (Bj + xl {j>n} ) C int E jt , i — l,...,k, E jt ^ E h for j { ^ j t } , (5.2) 

i.e., each of the sets (£>-,■ + xl{j >n }) is contained in a different cell of Tj. Using 
this notation we get 

lim P(T (1) n d = 0, T (1) n (d + x) = 0) 

>oo 

= V lim P(T« n d = 0, T« n (d + x) = 0, A,(x)) . 

Since the cells Si, 5 2 , . . . of Tj are finite with probability 1, we have 

lim P(T« n d = 0, T« n (C 2 + x) = 0, Aj(x)) = 

\x\— >oo 

if there are i < n and j > n with i, j £ Ji <E J. On the other hand, suppose 
that J = { Ji, . . . , Jfc} is a decomposition of {1 . . . , n + m} with Jj C {1, . . . , n} for 
z = 1, . . . , I and Ji C {n + 1, . . . , n + m} for i = Z + 1, . . . , k. Then we get that 

P(T« n d = 0, T« n (d + x) = 0, Aj(x)) 

= ¥(Aj(x), B, h n T« = 0, i = 1, . . . , Z, B Jt + x n T$ = 0, i = Z + 1, . . . , k) , 

where -Bj, = Uj^Bj and Tn y i, . . . ,Tn^ are independent copies of Tu which are 
independent of Tj. Thus we have 

¥(Aj(x), B Jt n r«. = 0, i = i, . . . , z, B Jt + x n t$ = 0, % = i + 1, . . . , k) 

= F(Aj(x))F(B Ji nT$ i = ®,i = l,...,l) P(^nT« =0,i = Z + l,...,fc). 
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Moreover, since Tj is mixing, we get 

lim ¥(Aj(x)) = P(Aj/(o)) P(Aj«(o)) , 

|a;|— >oo 

where J' = { Ji, . . . , J{\ and J" = {Ji+i, ■ ■ ■ , Jfc} are the decompositions of {1, ... , n} 
and {n + 1, ...,n + m}, respectively, induced by J, and Aji(o) resp. Ajh(o) are 
defined analogously to (15.21) . Summarizing the above considerations, we get 

lim P(T« D Ci = 0, T« n (C 2 + x) = 0, 

ir|^oo 

= F{A JI {o),Bj i nTff ti = Hi,i = l > ...,l) 

x P(Aj"(o), n Tff - = 0, i = I + 1, . . . , k) 
= P(T« D Ci = 0, Aj,(o)) P(T« n C 2 = 0, Aj„(o)) , 

which yields 

lim P(T« n d = 0, T« n (C 2 + ac) = 0) 

[as[— »oo 

= V lim P(T (1) n Ci = 0, T (1) n (C 2 + ar) = 0, Aj(x)) 

= ^ ^P(r( 1 'nc 1 = M J ,(op(TWnc 2 = Mj"W) 

j'eD' j"ez?" 
= P(T« n Ci = 0) P(T (1) n C 2 = 0) , 

where D', D" is the family of all decompositions of {1, ... , n} and {n + 1, . . . , n + m}, 
respectively. Thus, by Lemma I5TT1 the nested tessellation T is mixing. 



5.2 Integrability condition (13.31) 



The next result provides several classes of stationary tessellations such that the second 
moment of the circumference of their typical cell is finite, where i2(H) denotes the 
radius of the minimal ball containing the random convex polygon H. 

Lemma 5.3 If T is a PVT, PDT, PLT and STIT tessellation, respectively, then 
Ei? 2 (H*) < oo and, consequently, 

Ei^(9S*)<oo. (5.3) 
Moreover, (\5.3§ /io/ds if T is a a Tj/Tji- superposition /nesting such that 

max{ER 2 (E}),ER 2 (E* n )} < oo , (5.4) 
where 5} and S} 7 zs the typical cell ofTj and Tjj, respectively. 
Proof Note that 

Ev 2 (dE*) < 4tt 2 E j R 2 (S*) (5.5) 

holds for the typical cell H* of any stationary tessellation T. Furthermore, if T is a 
PDT, then it is well known that ER 2 (E*) < oo. This result goes back to [27], see 
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also Theorem 7.5 in [28] and Theorem 10.4.4 in [32]. Similarly, it is well known that 
ER 2 (E*) < oo holds if T is a PVT or PLT, see e.g. [7]. Since the interior of the typical 
cell of a STIT tessellation and a PLT have the same distribution ([30]). it is clear that 
ER 2 (E*) < oo also for STIT tessellations. If T = Xj/T/j is an iterated tessellation 
with cell intensity A^, then we can use Proposition 3.1 in [21] and Campbell's theorem 
in order to get 

Aj 



E / ^(d(E* n + x n 5*)) 2 { H| /+ ,nH^0} ^{dx) 



where A/, Xn and E}, E* n denote the cell intensities and the typical cells, respectively, 
of Tj and T/j. Note that we can assume that H| and E* n are independent random 
convex bodies. Since v\ {d{E* u + x H E})) < min{z/^(9S}), v\ (<95} 7 )} we get 

Ez/ 2 («9~*) < ^E^inl^ 2 ^)^^^)}^^©-*)) 

< f![^££ E ( mill{l ^ (SS j ))1 ^ (flE;j)} rnax^H*),^*,)} 

where in the latter inequality we used that 

is 2 {E* n ®E*) < vri? 2 (H*©H* 7 ) < 4vrmax{ J R 2 (S*), R 2 (E* n )} . 
Using (15.51) and the independence of E*j and S} 7 , this gives 

Evf(dE*) < E(min{iZ 2 (S}),iZ 2 (^ 7 )} max^S?),^,)} 



T 

"" > " E J R 2 (S 7 )E J R 2 (S 7/ ) < oo 



4vr 3 A / A // . 



provided that (15.41) holds. 

5.3 Asymptotic Weibull distribution of shortest path lengths 

In Sections [5 . II and I5T21 we showed that the assumptions of Theorem 13.21 are fulfilled for 
several classes of random tessellations T. Thus, we are now able to apply Theorem l3.2l 
to these tessellations. 

Corollary 5.4 Let Z ~ Wei(Xn,2) and let T be a PDT, PVT, PLT or STIT tes- 
sellation, or an iterated tessellation T = Tj/Tu such that condition < \5.^j is fulfilled, 
where T is either 

1. a superposition of two mixing tessellations Tj and Tu, or 

2. a nesting of a mixing initial tessellation Tj and any stationary component tes- 
sellation Tu. 
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Then C* —> £Z for some constant £ > 1 provided that 7 — > 00 and A^ — > suc/j £/ia£ 
A^7 = A. Furthermore, if T is a PLT or a Tj/Tjj -superposition/nesting, where Tj is 
a PLT, then £ = 1. If T is a PDT, then £ < 4/tt ss 1.27. 



Proof The first part of the assertion follows from Theorem 13.21 if the results of Lem- 
mas E2] and E3] as well as the comments immediately before Lemma I5T21 are taken into 
account. Now we consider the cases that T is a PLT, a Tj/T/j-superposition/nestrng 
with a PLT Tj, or a PDT. To begin with, let T be a PLT with intensity 1. Then, 
the edge set T 7 of the tessellation T 7 introduced in Section 12.51 is generated by a 
random sequence of lines L , L%, . . . , where Li,L 2 ,... form the edge set T 7 of the 
(stationary and isotropic) PLT T 7 and L is an isotropic line through the origin o, 
which is independent of T 7 . Thus we have 

-^ ) \^ n % r )) ^ i|*(T«nB(o,r)) + -. 



Using Theorem lA.il together with Lemma [A.21 this yields that the family of random 
variables {X 7j g, 7 > 0} with X 7 g = ^(T^T^nSfo/))/) is uniformly integrable 

since z/ 1 (T 7 1 ' 1 fl £>(o, r))/7 = Tcr 2 Ui(T^ fl B(o, ry))/u 2 (B(o, rj)) converges to r 2 n in 
L 1 due to the fact that the PLT T is mixing and, therefore, ergodic ([8], Theorem 
12. 2. IV). Furthermore, in Lemma EL~3l we showed that there is a £ > 1 such that 

p 

the Laplace transform of X lt ^ converges to 1, which implies that X 7 ,c — > (|20j, 
Theorem 5.3). Thus, applying Theorem IA.ll again, we get that 

limEX 7e = 0. (5.6) 

7—* 00 

However, if T 7 gets denser, there are lines which intersect the line through o close 
to o. Thus, all points on these lines have approximately the direct connections as 
shortest paths which can be used to show that £ = 1. Suppose that £ > 1 and let 
r > 2 > e > with £ > 1 + e. If the line Lj intersects the segment L e , where 
L ,e = Lq fl B(o,e/2), then for each y G Lj it holds that < c(y) — \y\ < £ since the 
path from y to o via the intersection point Li fl L 0i£ is not longer than \y\ + e. Thus, 
if \y\ > 2 

\c(y)-£\y\\ = \ c (y)- \v\ > e(l2/|-i)>e, 

which means that y G \T^ £ . Furthermore, if Li fl L , e 7^ 0, it is not difficult to 
see that ui(LiC\B(o, r)\B(o, 2)) > a for some constant a > 0. These two observations 
lead to 

^ = -^fW.nfl^r)) > -^i( |J {L.nsMVBM)}) 

> — : Li fl L , £ 7^ 0} 

7 

and, since : Li fl L , e 7^ 0} ~ Poi(2£y/Tc), 

liminf EX^ > lim - E#{L, : L» fl L , £ 7^ 0} = — > , 

7^00 ^ 7— >oo •y 7T 



17 



Scaling limits of shortest path lengths 



which is a contradiction to (I5.6p . Thus, £ = 1 holds. If the tessellation T = Tj/Tjj is 
a superposition/nesting such that Tj is a PLT, then 



" H B(o, r)) > I {oef a) } - ^(T«\T« £ n S( 0> r)) 



where Tj' denotes the part of T) which corresponds to Tj. Since Tj is assumed to 
be a PLT, the same arguments as above can be applied to show that £ = 1. Finally, 
let T be a PDT and let N(y) denote that node of T which is closest to y G M. 2 . It 
has been shown in [4] that for any t > and y G dB(o, 1), there is a path P{ty) from 
N(o) to N(ty) on T"W with length c(P(ty)) such that almost surely 

lim C -^l = A - . (5.7) 

t >oo t IX 

Consider the stationary point process n L of intersection points {Xj}, where 
L = {sy : s G R} and • ■ ■ < X_ x < X < < X x < ■ ■ ■ , and denote by c(X i5 Xj) 
the shortest path length from Xj to on T^ 1 ' . Furthermore, consider the stationary 
marked point process {(Xj, c(Xj, iV(Xj))} and denote its typical mark by c^. For 
each i > we then have 

c(X ,X) < c(N(o),N(X t )) | c(X ,X(o)) | c(X,X(X)) 
|Xj — Xo| |Xj — Xo| |Xj — Xo| |Xj — Xo| 

< c(P(X)) | c(X ,X(o)) | c(X,X(X)) 
|Xj| |Xj| |Xj] 

Clearly, the second summand of the latter expression tends to as i — ► oo. The same 
is true for the third summand, because {(Xj, c(Xj, X(Xj))} is ergodic and Ec^ < oo. 
Thus, by (15.71) . we get that 

r c(X ,X) 4 

hmsup— — — < — . (5.8) 

i->oo |Xj — Xo| 7T 

On the other hand, we have P(limj^ 00 c(X , Xj)/|Xj — X | = £) = 1 if and only if 
P(lim^ oo c(X *,X;)/|X* - X*| = f) = 1, where {X*} is the Palm version of {X}. 
Now, using (EHD and (EH), it follows that £ < 4/vr. 



5.4 Some extensions 

Note that the setting of Theorem 13.21 can be generalized in different ways. For exam- 
ple, the statement of this theorem remains valid if instead of C* the typical subscriber 
line length S* is considered, where S* is the shortest path length from the origin to the 
nearest point X#,o of X#, which is defined as the sum of the distance from the origin 
to the nearest point of the edge set and the shortest path length on from 
this point to Xh,o ( [1 3j ) . Note that in this case the auxiliary results corresponding 
to Lemmas 14.31 and 14.41 can be proved basically in the same way. 

Furthermore, in the proof of Theorem 13.21 it is not necessary to assume that T is a 
random tessellation, but it is possible to consider an arbitrary stationary and isotropic 
segment process in IR d which is mixing and such that there is only one single cluster 
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(a) k — ► (where \i = 1) (b) k — * oo (where j\e = 1) 

Figure 3: Densities of C* if T is a PLT (together with corresponding limit distribu- 
tions) 



with probability 1. This means in particular that Theorem 13.21 can be extended to 
random geometric graphs. 

Another kind of extensions can be obtained by relaxing the assumption that X^n 
is connected to the nearest point of Xh, i.e., T# is a Voronoi tessellation. For instance, 
Xl, u can be connected to its fc-th nearest neighbour of Xh for any k > 1. Then, in 
Theorem 13.21 we only have to replace Z by the distance from the origin to the k-th 
nearest point of a Poisson process which is distributed according to a generalized 
Gamma distribution ([151 SI])- Further possible extensions include that T# is a 
certain Cox-Laguerre tessellation Q23]) or an aggregated tessellation ([21 134]). 



6 Conclusion and Outlook 

We consider the typical shortest path length C* of stochastic network models with 
two hierarchy levels, where the locations of network components are modelled by Cox 
processes on the edges of random tessellations. It is shown that the distribution of 
C* converges to known limit distributions for extreme cases of the model parameters, 
i.e., if a certain scaling factor k tends to zero or infinity. 

The results of the present paper have applications in the analysis of telecommu- 
nication access networks since the distribution of C* is closely related to cost and 
risk analysis of such networks ([14]). Using the fitting techniques introduced in [1 2j , 
an optimal tessellation model can be chosen for a given set of road data. Moreover, 
the scaling factor k can be estimated. Then, on the one hand, for small values of k 
the limit distribution of C* is directly available and it does not depend on the type 
of the optimal tessellation model. On the other hand, for large values of k the limit 
distribution of C* and an upper bound for this distribution is directly available if the 
optimal model is PLT or PLT-superposition/nesting and PDT, respectively. 

In order to get an idea how small or large the scaling factor k should be (to 
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replace the distribution of C* by the corresponding limit distribution) and how to 
calculate the constant £ appearing in the limit distribution for C* as k —* oo, the 
density of C* can be estimated by Monte Carlo simulation of the typical serving zone 
(|36j). This can be done for PVT, PLT and PDT as well as for superpositions and 
nestings built from these basic tessellation models, using simulation algorithms of 
the typical serving zone introduced in [9l dU ETJ GS|- I n Figures [3] and 2] estimated 
densities for different values of k are shown together with the corresponding limit 
distributions if the tessellation model chosen for the underlying road system is a 
PLT and PVT, respectively. As can be seen in Figure 0] (b), the density of the 
Wei(\n/1.14:5 2 , 2)-distribution approximates the density of C* very well for T being 
a PVT and k > 1000. This suggests that in this case the constant £ appearing in 
Theorem 13.21 and Corollary 15.41 respectively, is approximately 1.145. 

Furthermore, the limiting distributions derived in the present paper can be used 
to choose parametric densities which can be fitted to the estimated density of C* for a 
large range of k. Parametric families which include both exponential distributions and 
Weibull distributions turned out to be good choices, see [14]. In Figure [5] estimated 
densities for different values of k are shown together with fitted truncated Weibull 
distributions. Note that these truncated Weibull distributions have two parameters 
and there is a quite good fit for both tessellation models considered in Figure [5] and 
for a large range of values of k. 



A Some mathematical background 

In the proof of Lemma l4~3l given below we make use of some well-known results from 
measure theory, the theory of subadditive processes, and geometric measure theory 
which are briefly summarized. We start with the definition of convergence in measure 
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(a) PLT (b) PVT 

Figure 5: Densities of C* (together with fitted parametric densities) 



and uniform integrability which can be used to characterize //-convergence. A family 
of measurable functions {/ 7 , 7 > 1} defined on a measurable space (fl,A,fi) and 
taking values in R converges locally in /z-measure to a measurable function / : O — > R 
if 

lim/x({|/ 7 -/|>£}nA)=0 (A.l) 

for all e > and A 6 A with fi(A) < 00, where \i is assumed to be a a-finite measure. 
If fj, is a probability measure such that (IA.1I) holds for each e > and A = Q, then 
one says that / 7 converges in probability to /. Furthermore, if for each e > there 
is a /i-integrable function g such that 

/ \f 7 (u})\fji(du}) <e for all 7 > 1 , (A.2) 

h\fi\>9] 

then the family {/ 7 , 7 > 1} is said to be uniformly /x-integrable. With the above defi- 
nitions it is possible to characterize the //-convergence as follows; see Theorem 2.12.4 
in ©. 

Theorem A.l A sequence of fi-integrable functions fx, $2, ■ ■ ■ '■ — > R converges in 
L 1 to a n-integrable function / : O -» E i/ and on/y «/ (i) /„ converges locally in 
jjL-measure to f and (ii) {/ n } is uniformly [i-integrable. 

We still mention an elementary but useful result which immediately follows from the 
definition of uniform integrability. 

Lemma A.2 Let {/ 7 ,7 > 1} and {g 7 , 7 > 1} 6e two families of measurable functions 
on (ft, A, //) which satisfy that |/ 7 | < |p 7 | /or a// 7 > 1. T/«en {/ 7 ,7 > 1} is uniformly 
fi-integrable i/{g 7 ,7 > 1} is uniformly fi-integrable. 

Another useful tool is the notion of subadditivity. Let Y be a family of real- valued 
random variables Y = {1^,0 < i < j} defined on some probability space (Q,A,F). 
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Note that Y can be seen as a random element of some measurable space (S, B(S)) of 
double-indexed sequences, where B(S) is the Borel-a-algebra of S. Then Y is called 
a subadditive process if 

1. Yik < Yij + Yj k for all i < j < k, 

2. Y = {y -}iY' = {y (+lj+1 } I 

3. EF | < oo, where = max{0, loi}- 

The following result is due to Kingman ([21], Theorem 1). It is called the subadditive 
ergodic theorem; see also Theorem 10.22 in |20j . 

Theorem A. 3 Let Y be a subadditive process. Then the limit 

C = lim -Y 0] (A.3) 

exists and is finite with probability one and E£ = mfj^ILYoj/j. IfE£ > — oo, then 
the convergence in (\A.5ft) also holds in the L x -norm. Moreover, let 1$ C B(S) be 
the a-algebra of subsets of S which are invariant under the shift Y i— > Y', where 
YL = Yi + xj + i, and let 1 = Y^Xs C A be the corresponding sub-a- algebra of events. 
Then, 

C = lim - E(Y 07 - I J) . (A.4) 

Note that a subadditive process Y is called ergodic if P(Y G A) = or P(Y G 
A) — 1 for each A G T s . Thus, in the ergodic case, the limit ( considered in (|A.3I) 
and (IA.4I) . respectively, is almost surely constant. 

Finally, we use a decomposition of the Hausdorff measure V\ which is a special 
case of the generalized Blaschke-Petkantschin formula ([IB], Proposition 5.4). 

Theorem A.4 Let C C M 2 be a differentiable curve and assume that 

vi({x G C : Tan[C,x] = span{x}}) = , (A. 5) 

where Tan[C, x] is the tangent at x to C and span{x} = {cx : c G M} is the line 
which goes through the origin o G M 2 and the point x G C. Then, for any measurable 
g : C — > [0, oo) it holds that 

I g{x)u x {dx) = / V -J^L^)^, (A.6) 

where L\ is the half line of direction $ G [0, 2ir) emanating from o and CKj is the angle 
between Tan[C,Xj\ and span{xi}. 
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B Proof of Lemma 14.3 



With the help of Theorems lA.ll - IA.4l stated above we are now able to prove Lemma l4~3l 
Obviously, lim sup 7 _ too E exp v x {f\ x) \T^ e nB(o, r)) ) < 1. Thus it is sufficient 
to show that 

A 

7— >oo \ 7 



lim inf E exp ( - - v x (T« \ T^l e nB(o, r)) 

7— >oo \ ry v ' 7>?i 6 1 



> 1 . 



(B.l) 



Proof of ( IB. If ). First recall that we can identify with the Palm version A* (1) 

of the stationary random measure A (1) given by h T m(B) = v\(B PI Tj) for B e B 2 
since A (1) is the random driving measure of the Cox process X L , see [33], p. 156. 
Then, using the abbreviation 

A 

1 



h(rW) = exp(-- z/!(r« \ rg n S(o,r))) , 



where r^* ={m6 : |c(w) — £\u\\ < e} and c(u) denotes the length of the shortest 
path from u to the origin along the edge set of a tessellation r with o 6 r^, we get 
from the Campbell theorem for stationary random measures ([8], Proposition 13. 2. V) 
that 

where we used the substitution z = 72 in the last expression bearing in mind that 
(1/7)T"W = Ty . Furthermore, we put T^} iZ = {y G : \c(y, z/j) — £\y — z/j\\ < 
e}, where c(y,z/j) denotes the length of the shortest path from y to 21/7 along the 
edges of the considered graph. Then, for each 7 > 1, we get that 



E/>(T«) =-e( I exp(--u l {T^\T^lnB(zh,r)))u 1 (dz) 

> - EU(T^ n B(o, 1)) inf exp f - - ^(T«\T« 2 n B(z/ 7 ,r)) 



= - EL(r« n B(o, 1)) exp ( - sup - Vl (T^\T^ z fl B{z/y,r)) 

> i E(Wt« n 5(o, 1)) exp ( - sup - ^(rW\T« 2 n 5(o, r + 1)))) . 
^ v v «eT(i)ns(o,i) 7 /y 

Now, in order to prove (IB.ll) . it is sufficient to show that 

^ for 7 -> 00, (B.2) 

where X 7j ,e = sup zeT (i) nB(0;1 ) ^ 1/1 ^T^^T^z fl B(o,r + 1)Y To see this, note first 
that (IB.2I) implies that X 7) g converges in probability to 0. Thus, the random variable 
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Y 7i5 = exp(-AX 7i5 )z/ 1 (T( 1 ) n B(o,l)) converges in probability to V\ (T^ n B(o, 1)) 
if (HJ) holds. Moreover, < v x {TW n S(o, 1)) for all 7 > 1 and Eu^T^ n 
_B(o, 1)) = 7r < 00, which means that {Y 7 ^, 7 > 1} is uniformly integrable. Hence, 
Theorem I A. II yields that Y 7 £ converges in L 1 to u x (T^ f)B(o, 1)) and, in particular, 
lim^oo 1/ttEY 7i £ = l/nEu^T^ n B(o, 1)) = 1 if {EH) holds. Thus, jEjJ follows if 



we can show that (B.2) is true. 



Proof of (IB. 2D . Since X^g > it suffices to show that E 7) £ — > 0. Furthermore, 
note that with probability 1 the segments of the segment system Tj R B(o,r + 1) 
fulfill the conditions of Theorem IA.41 since none of these segments „points" to the 
origin because was assumed to be stationary. Thus, using Theorem IA.4I we get 
that 

EX 7iC =e( sup -f %^ ) (\c{y : -)-i\v--\\)v 1 {dy)) 
K zeTWnB(o,i) 7 JT^nB(o,r+i) V 1 7 7 7 J 



GT( 1 )nB(o,i) 7 ir^nB(o,r+i) 

1 f* 27T I ~W~ I 

e(- sup / v (| c ( Xj ,-)-eix--i|)d$) 



|Xi|<r+l 



<—e([ T sup V -J-l^dc^,-)-^--! 

7 Vo zG T(i)nBfo,l) m sma i V 7 7 



/o zeT( 1 )n-B(o,i) XjeT a) ni +. 

|Xi|<r+l 

= 2 ! r(r ± l) E / ^ -J-l^flc^^-^-il 

7 \ e r(DnB( ,i) V smai 7 7 

|X t |<r+l 

= 27r(r + l)E^ 7 (T( 1 )), 

where in the last but one line we used Fubini's theorem and the isotropy of T 7 , 
denoting by L + = Lq the half line with direction $ = 0, and in the last expression 
we used the abbreviation 

S ,.(T< 1 ») = i sup V -^—l M (\c(X i ,-)-(\X i --\\). (B.3) 

|Xj|<r+l 

Since the point process T^nl is stationary with intensity 2/n (|32j. Theorem 4.5.3), 
where we identify R with the x-axis, we can apply the inversion formula for Palm 
distributions of stationary point processes on R; see Proposition 11.3 (iii) in |20j. 
Thus, if T^* denotes the Palm version of with respect to the point process 
T (1 ) n R, we get that 

E^ 7 (T«) = -E / l [0 , xl] (x)g 7 (T^*-x)dx), 
n v Jo 7 

where the points of {X*} = T^* fl R are numbered in ascending order such that 
. . . < X* x < Xq = < Xj* < X2 < ■ ■ .. Hence, in order to prove (|R2i) it suffices to 
show that 



(•00 

lim E ' 

7^00 



( jf I [0l x ?] (x) <7 7 {T {1) * - x) dx) = , (B.4) 
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where the function g 7 : T — > [0, oo) is given in (IB.3I) . The proof of (IB.4I) is subdivided 
into two main steps. First, we show that 

lim gJx,T {1 >) = (B.5) 

7— >oo 

almost everywhere with respect to the product measure V\ <8> P*, where we used 
the abbreviating notation # 7 (x, T"W*) = 1[ 0>x *](x) g 1 {T^ v> * — x) and P* denotes the 
distribution of T^* . Then, we show that {g 7 , 7 > 0} is uniformly (z/i®P*)-integrable. 
By means of Theorem IA.ll this implies that (IB.4I) holds. 

Proof of rtBJjj ). Note that for each x G [0, X*] we get 

< - sup V -^—%J\c{X h -) - %\Xi - -I |) 

\Xi\<r+l 

= - sup ^ -^i^f-l^x^^-eix;-^!) 

7 z( zTW*riB!r 11 , sinaj V7 1 v 

X*£B(x,(r+l'h) 

<- Y\ — sup l [£t00) (-\c(X*,z)-Z\X*-z\ 
7 f— ' sm aj zP TW*nB(T 11 ^7 

X*eB(x,(r+l)j) 

Thus, 

|3ff]<(r+a) 7 

= - T^— I[ E ,oo)(- SUp | C (X*,z)-£|X*-2||), 

7 x , gT ^ ni+ sma « V 7 «TWnBM 7 

|X;|<(r+a) 7 

where a = 1 + X*. Furthermore, we have 

- sup \c{Xlz)-t\X*-z\\ < - (c(o,X*)-Z\X*\) 
7 zeTm*nB{o,a) 7 

+ — ( sup c(z,o)+C,a 

7 \z€TW*nB(o,a) 

since c(X*, o) - c(o, 2) < c(X*, 2) < c(X*, o) + c(o, 2) and - £a < - z\ < 

£|X*|+£afor alii > 1 and z G T«*n J B( o, a). Clearly, the second term of this upper 
bound tends to zero P*-almost surely as 7 — > 00. Thus in order to show that (jB.ol) 
holds, it suffices to prove that P*-almost surely 

I(c(o,X;)-W)g(-|, I) (B.6) 
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for all sufficiently large i > 1 such that X* < (r + 0)7. 

Proof of (IB. 61) . Note that X = {|X* - X*\,i,j > 1, i < j} is an additive process, 
because |X* - X*| = \X* - X*\ + \X* - X*| for 1 < j < k. Since T«* n R is 

cycle-stationary (see e.g. [35]), we have that {|X* — X*|} — {\X* +1 — X* +1 |}, where 
< KXl < 00. Thus, by Theorem IA.3I we get that the finite limit lim^oo X* Ji = Cx 
exists P*-almost surely. Furthermore, consider the family Y = {Yij,i, j > l,i < j} 
of non-negative random variables with = c(X*,X*), where c(X*,X*) denotes the 
shortest path length from X* to X* on T"W*. Then, it is easy to see that < Yij+Yjk 

for i < j < k. By the cycle-stationarity of T^* fll, we have that {Yij} — {Y i+ ij + i}, 
where EF i = Ec(Xq,Xi) < 00 holds by condition (13.31) : see the next paragraph 
below. Thus Y is a subadditive process and we can again apply Theorem IA.3I to get 
that the finite limit linx^oo c(Xq, X*)/j = Cy exists P*-almost surely. Since X and Y 
are ergodic (see the paragraphs below), the limits Cx and Cy are constant. Noticing 
that < EX* = Cx < Cy < 00, this gives that 

h ^ =Il i« ={ , (B.7) 

where £ = Cy/Cx e [l,oo). Now let e > such that e(r + a) < e/2. Then ijETrj) 
implies that with probability 1 

for all i sufficiently large and, therefore, 

x - {c{o,x;)-ix;) e (- 1 , I 

if i is sufficiently large and X//7 < r + a. 

Proof of Ec(Xg , X*) < 00. Consider the stationary marked point process {(X n , 5+)}, 
where {X n } = n R is the point process of intersection points of the edge set 
with the line R, and S+ the cell of T on the right of X n . Let A + denote the intensity 
of the marked point process {(X n ,S+)}, and S + * its typical mark. Then, by the 
definition of the Palm mark distribution (see e.g. Section [2T2T) . we get that 

Ec(X*,X*) < Ei^SE**) = ^E Y, U ^ 9E ^ 

XjeTWnlo.i) 

= E S I {9+H l n[O,l)^0} Z/ l(^) > 

where <9 + H denotes that part of the boundary of 5 with outer unit normal vector in 
[ir/2, 37r/2). Thus, applying Campbell's theorem to the latter expression, we have 



Ec(X *,X*) < ^Eu^dET) [ 1 



{3+S*+a;n[O,l)^0} 



u 2 (dx) 



^E i / 1 (aH*)^ 2 ([0,l)©a+H*) 
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where A T = 1/Ez/ 2 (H*). Since u 2 ([0, 1)©<9 + H*) < avi(&E*) for some constant a < oo, 
this implies that Ec(Xg,X*) < (a\ T /X + ) Eisft&E*). Thus, the assertion is shown. 

Ergodicity. We only prove that X is ergodic, because the ergodicity of Y can be 
shown in the same way. Recall that by X s C B{S) we denote the cx-algebra of those 
subsets of the space S of double-indexed sequences, which are invariant under the 
shift {\X* - X*\} i — ► {\X* +1 - X* +1 \}. Furthermore, note that X = h(T^*) for 
some measurable function h : T — > S, where for any tessellation r in M. 2 and A G Is, 
we have h(r^) E A if and only if h(r^ — x) E A for all x E [0, oo). Thus, from the 
definition of the Palm distribution of the stationary point process {Xi} = T"W n E 
with intensity 2/tt, we get for any AeX s that 



^E J2 UK^-Xtf 



X i GT( 1 )nB(o,l)nL+ 

71 



{XEA) = P(/i(T (1) *) E A) 

IT 

2 

2 E(l A (h{T^)) E T« n B(o, 1) n L+}) 

= I E(l h -i (A) (T«) G rw n B(o, i) n L+}) . 

On the other hand, since T\ is mixing and h~ x {A) = h^{A) + x for any A G Xs and 
x G L + , we have 

P(T (1) G /T 1 ^)) = lim P(T (1) G /T 1 ^), T (1) - z G /T 1 ^)) 

= P(T (1) G /^(A)) 2 , 

which implies that PfT^G /i _1 (,4)) = or P(tMg fc -1 ^ 4 )) = 1- Th us, altogether, 
we have 

p(XgA) = P(T (1) er 1 (A))^E#{l i eT (1) nB(o 1 i)nL + } 

= P(T (1) G 

and, consequently, P(X G A) = or P(X E A) = 1 for any A E X s , which means 
that X is ergodic. 

Uniform integrability. Finally, we show that the family {g 7 , 7 > 0} considered 
in (IB.5I) is uniformly (v\ ® P*)-integrable. From the ergodic theorem for stationary 



marked point processes ([8], Theorem 12. 2. IV), we get that 

lim — = (r + 1) lim — 

7^007 sina,- 7-^00 (r + 1)7 sina,- 

|Xi|<(r+l)7 \Xi\<(r+l)-y 

= (r + 1) E(sina*)" 1 

almost surely and in L 1 since the point process fl M marked with the intersection 
angles is ergodic, which can be shown in the same way as the ergodicity of X. Here 
a* denotes the typical intersection angle which is distributed according to the density 
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f a ,(a) = sin(a)/2 for < a < vr, see e.g. [33], p. 288. This yields E(sina*)- 1 
7r/2 < oo. Thus 







lim E 

7^00 



- V — (r + nEfsina*)- 1 

-v ^-^ sin n/- ' 



7 * — ' sm a 

X 4 GT( 1 )nL+ 
|X i |<(r+l) 7 



lim - E / I [0jXI] 



X 



1 

7 



Xi&iT^* -x)r\L+ : 
l*i|<(r+l)7 



sm ola 



[r + 1) E(sina 



dx . 



where in the last equality we used the inversion formula for Palm distributions of 
stationary marked point processes on R; see Proposition 11.3 (iii) in [2D]- In other 
words, we showed that 



2 



lo,xi]( x ) ~ Y] —— 

1 7 ' sm«i 

X i e(T<. 1 >-x)nL+: 
|Xi|<(r+l) 7 



+ l)I [0 ^ ?] (a;) E(sina*) 



in L x {y\ ®P*) as 7 — > oo. This means in particular that the family {/i 7 , 7 > 0} with 
h,(x,TW*) = l [0 , xl] (x) I- ' 



7 < — ' sm a; 

X l £(TW*-x)llL+: 
|Xi|<(r+l)7 



is uniformly {y\ <g> P*)-integrable; see Theorem IA.11 Furthermore, we have that 
l [0 , xl] (x)g 7 (TW* -x)< %* f] (s) ^ ^ ' 



7 



X i e(T( 1 )*-a;)nL+: 
l^i|<(r+l)7 



sm a,- 



Thus Lemma IA.2I yields that the family {g 7 , 7 > 0} considered in (1B.5I) is uniformly 
(v\ <8> P*)-integrable. 
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